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The neural crest is an embryonic stem cell population that gives rise to a multitude of derivatives. In particular, the cranial 
neural crest (CNC) is unique in its ability to contribute to both facial skeleton and peripheral ganglia. To gain further 
insight into the molecular underpinnings that distinguish the CNC from other embryonic tissues, we have utilized a CNC- 
specific enhancer as a tool to isolate a pure, region-specific NC subpopulation for transcriptional profiling. The resulting 
data set reveals previously unknown transcription factors and signaling pathways that may influence the CNCs ability to 
migrate and /or differentiate into unique derivatives. To elaborate on the CNC gene regulatory network, we evaluated the 
effects of knocking down known neural plate border genes and early neural crest specifier genes on selected neural crest- 
enriched transcripts. The results suggest that ETSI and SOX9 may act as pan-neural crest regulators of the migratory CNC 
Taken together, our analysis provides unprecedented characterization of the migratory CNC transcriptome and identifies 
new links in the gene regulatory network responsible for development of this critical cell population. 



[Supplemental material is available for this article.] 

The neural crest is a multipotent cell population that contributes 
to a wide variety of derivatives, including sensory and autonomic 
ganglia of the peripheral nervous system, cartilage and bone of the 
face, and pigmentation of the skin (Le Douarin and Kalcheim 
1999). Acquisition of the neural crest played a critical role in ver- 
tebrate evolution, facilitating evolution of the jaw and active 
predation that fostered diversification of vertebrates (Gans and 
Northcutt 1983). Due to the extensive contributions of the 
neural crest to vertebrate embryos (Hall 2000), defects in its de- 
velopment contribute to a large percentage of congenital birth 
defects (Trainor 2010). 

Neural crest cells are induced during gastrulation (Basch et al. 
2006), but first become apparent as "premigratory" cells in the 
dorsal neural tube by expression of a suite of neural crest specifier 
genes like FOXD3, SNAI2, and SOX9 (Sauka-Spengler and Bronner- 
Fraser 2008). They subsequently undergo an epithelial to mesen- 
chymal transition (EMT) to leave the neural tube, migrate long 
distances, and differentiate into diverse cell types. To understand 
the progression of their development, we have proposed a gene 
regulatory network (GRN) that underlies this complex process 
(Meulemans and Bronner-Fraser 2004), comprised of nested sub- 
networks, each controlling distinct events from specification to 
emigration, and ultimately to migration and differentiation 
(Betancur et al. 2010a). However, this network is far from com- 
plete: A relatively small number of transcriptional regulators have 
been identified/characterized, and little is known about effector 
genes present in migrating neural crest cells (Simoes-Costa and 
Bronner 2013). Thus, elaboration of the neural crest GRN is cru- 
cial to clarify how distinct differentiation programs are activated 
in multipotential precursors to drive commitment toward specific 
fates. 
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Although all neural crest cells share common molecular and 
behavioral traits, they can be subdivided operationally according 
to their axial level of origin, as cranial, vagal, trunk, and sacral (Le 
Douarin and Kalcheim 1999). These subpopulations differ in mi- 
gratory pathways and types of derivatives formed; e.g., cranial 
neural crest (CNC) has the ability to give rise to cartilage and 
bone, whereas trunk neural crest does not, even after trans- 
plantation in the head (Nakamura and Ayer-le Lievre 1982; Le 
Douarin et al. 2004). Recent evidence points to distinct regu- 
latory programs controlling expression of neural crest specifier 
genes SOX 10 and FOXD3 in these two subpopulations (Betancur 
et al. 2010b; Simoes-Costa et al. 2012), raising the intriguing 
possibility that a molecular framework may underlie such 
differences. 

Here, we take advantage of a cranial-specific neural crest en- 
hancer as a tool to isolate pure populations of migrating CNC cells 
from an in vivo context, allowing us to obtain a genome-wide 
representation of the active transcriptome of this important cell 
population. Because migrating CNCs intermingle with other 
mesenchymal cell types, obtaining pure populations has been 
particularly difficult. Moreover, no available reagents distinguish 
CNCs from subpopulations at other axial levels. Recently, how- 
ever, we have isolated NC enhancers that mediate gene expression 
at unique stages and axial levels (Betancur et al. 2010a), including 
SoxlOE2, which mediates reporter expression selectively in all 
migrating CNC. This tool enables us to isolate these cells and ob- 
tain genome-wide information regarding their active tran- 
scriptome using RNA-seq. Here, we utilize this approach to char- 
acterize the suite of genes enriched in migratory CNCs. The results 
provide unprecedented characterization of their transcriptome 
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and allow identification of new links in the GRN responsible for 
development of this critical cell population. 

Results 

Isolation of pure cranial neural crest population 

Expression of the specifier gene SOX10 initiates as NC cells com- 
mence emigration from the neural tube at all axial levels (Kim et al. 
2003; McKeown et al. 2005). Previously, we identified avian NC 
enhancers that drive reporter expression in a manner that re- 
capitulates endogenous SOX10 expression (Fig. 1A; Betancur et al. 
2010b). In particular, SoxlOE2 (Fig. 1B,C,C) mediates expression 
in the CNC but not in the trunk subpopulations. 

To confirm specificity, we compared enhancer activity with 
endogenous SOX10 expression. Antibody staining in embryos 
electroporated with Soxl 0£2:eGFP shows that the enhancer is only 
active in SOX10+ neural crest (Fig. 1D-F). Although the majority of 
ectodermal cells are efficiently coelectroporated with SoxlOE2 en- 
hancer-driven (eGFP) and ubiquitous promoter-driven (mRFP) 
constructs, as highlighted by H2B-RFP expression (Fig. 1G), only 
SOX10+ cells exhibit eGFP expression (Fig. 1H-J). To examine en- 
hancer activity during neural crest migration, we performed time- 
lapse imaging. Supplemental Movie 1 highlights the remarkable 
movement of CNC cells as they migrate throughout the embryonic 
head (Fig. 1K-M). Thus, SoxlOE2 enhancer is a useful tool for in vivo 
labeling and identifying migrating CNCs. 

In order to obtain a transcriptional profile of the CNC, 
SoxlOE2:eG¥P was electroporated into gastrula-stage chick em- 
bryos subsequently incubated to stages of active CNC migration 
(HH10) (Fig. 2A). After dissection and dissociation of embryos, 
a pure population of migrating CNC cells was isolated by FACS 
with 7-AAD exclusion to eliminate unhealthy or dying cells. Ap- 
proximately 5%-10% of the viable dissociated cranial cells 
expressed eGFP, resulting in about 500-1000 labeled CNC cells per 
electroporated embryo. After sorting, cells were tested for enrich- 
ment of known neural crest markers such as SOX1 0 and SNAI2 (Fig. 
2C) using qPCR. Isolated GFP+ sorted CNC cells had strong en- 
richment for both genes relative to GFP- cells. Subsequently, RNA 
of the purified CNC population was extracted and reverse tran- 
scribed. cDNA was amplified using a linear amplification system 
and used for sequencing library building. 

Identification of genes enhanced in migrating cranial 
neural crest 

The CNC RNA-seq results were compared with those of GFP- cells 
sorted from whole embryos of the same stage to identify transcripts 
enriched in migrating CNC. The results revealed numerous mole- 
cules strongly enriched in migrating crest compared with the 
whole embryo. We identified 1301 significantly up-regulated 
genes (Fig. 3A; Supplemental Table 1) that satisfied criteria we 
established for the CNC data set (see Methods). Enhanced genes 
included phosphatases, transporters, peptidases, transmembrane 
receptors, kinases, ion channels, G-coupled receptors, transmembrane 
receptors and transcriptional regulators. Figure 3B shows the rel- 
ative abundance of transcripts grouped by molecular function. 
Enzymes accounted for 41% of up-regulated transcripts, followed 
by transcriptional regulators (17%) and transporter molecules 
(11%). The data set also included a large number of uncharac- 
terized proteins (83 genes) (see Supplemental Table 1), perhaps 
reflecting limitations of the current chick genome assembly. In 



addition, 257 genes were depleted in the migrating CNC pop- 
ulation relative to the whole embryo (Supplemental Table 2). 

As proof of principle regarding data set quality, we looked for 
known transcription factors, including bona fide neural crest 
markers such as TFAP2A, PAX7, SOX8 and SOX 10, FOXD3, among 
others. As expected, all were present in our data set (Fig. 3C). 
SOX 10 is one of the genes with highest enrichment relative to 
GFP- cells (Fig. 3D). ETS1, a key regulator of CNC formation 
(Theveneau et al. 2007; Betancur et al. 2010b; Simoes-Costa et al. 
2012), is significantly enriched in the data set although it is 
strongly expressed in other embryonic tissues such as endothelial 
cells. 

Figure 3D shows genes with high enhancement or depletion 
in the migratory CNC. In addition to transcription factors, genes 
known to play important roles in CNC formation are also present. 
For example, Collagens 2 and 9, crucial for chondrocyte de- 
velopment, are enriched (Fig. 3D), as are nuclear receptors such as 
RXRG and other transcriptional regulators such as ALX1. Genes 
visibly down-regulated or absent from the CNC include a number 
of HOX genes, which is expected since the CNC lacks most HOX 
genes, GATA factors 4,5,6 as well as TBX5, 20 and 15 (Fig. 3C; 
Supplemental Table 2). Importantly, ZIC1, important for trunk NC 
specification and a trunk-specific input in FOXD3 (Simoes-Costa 
et al. 2012), is also absent from the CNC (Fig. 3C). 

Biological pathway analysis of cranial neural crest 
transcriptome 

To gain an overview of the molecular processes operational in the 
CNC transcriptome, we performed a biological pathway analysis of 
the data set. Confirming its utility, the results reveal enrichment in 
genes involved in pathways known to be important for neural crest 
migration, such as integrin (Fig. 4A) and CXCR4 signaling (Fig. 4B). 
Revealing predictive value, the analysis also suggests molecular 
processes that have yet to be investigated in depth during neural 
crest formation. For example, a large number of molecules related 
to RNA processing are up-regulated in the CNC (Fig. 4C). This 
suggests that RNA-binding proteins and differential splicing 
mechanisms may be important in neural crest development. 
Similarly, a number of molecules involved in post-translational 
modification, such as the ubiquitination pathway, are present in 
the data set. This result is consistent with recent findings in- 
dicating the importance of molecules involved in post-trans- 
lational modification in regulating aspects of neural crest de- 
velopment (Lee et al. 2012). The analysis also points to a potential 
role for pathways that have yet to be investigated in neural crest 
development, such as the IGF1 and JAK/STAT pathways (Fig. 4E). 

A number of genes in the CNC transcriptome data set are 
disease-related. Our analysis revealed association of CNC genes 
with diseases such as osteochondrodysplasia (Fig. 4D), consistent 
with the importance of CNC for the formation of the skeletal el- 
ements of the face. Furthermore, the data set contains genes linked 
to malignant diseases, such as neuroblastoma, glioma, and other 
types of neural crest derived cancers (Fig. 4E). This is expected due 
to the similar behavior of neural crest cells and metastatic cancers 
and the fact that several neural crest genes such as MYC and CXCR4 
are routinely used as biomarkers for tumor diagnosis. 

Novel transcriptional regulators for the cranial neural crest 

To validate the CNC transcriptome data set, we performed a sec- 
ondary in situ hybridization screen (Fig. 5). The results reveal many 
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Figure 1. SOX7 0c/s-regulatory element drives specific reporter expression in the migratory cranial neural crest. (A) Neural crest specifier gene, SOX10, is 
expressed by all migrating neural crest cells. (8) A chicken embryo electroporated with Soxl 0E2 enhancer has GFP expression in migratory CNC at HH1 0. 
(C,C) Transverse sections of HH10 chicken embryo electroporated with 5ox70£2:eCFP plasmid. GFP+ cells are also immunoreactive for the HNK-1 
epitope. (D-F) Colocalization of GFP+ cells and endogenous SOX10 confirms specificity of Sox 1 0E2 enhancer in the CNC. (G-J) Confocal imaging confirms 
that of all electroporated cells (red), only the ones with SOX10 positive nuclei (blue) express GFP. (K-M) Time-lapse imaging of a chicken embryo 
electroporated with 5ox70£2:eGFP demonstrates the extensive movements of cranial neural crest in the embryonic head (see Supplemental Movie 1). 
(NC) Neural crest; (R3) rhombomere 3; (OV) otic vesicle. 



interesting candidate genes that may play a role in the process of 
CNC migration and/or differentiation. Eighty-five genes were 
cloned and their expression patterns verified by in situ hybrid- 
ization. Approximately 90% of genes had strong expression in the 
migrating CNC as expected (Supplemental Table 3). We also tested 



enrichment of a number of genes by RT-qPCR in GFP+ cells; results 
are consistent with the RNA-seq data (Fig. 5P). 

Our in situ hybridization screen shows specific CNC gene 
expression for molecules strongly up-regulated in our data set 
(Fig. 5). Genes such as transcription factor EBF1, adaptor protein 
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Figure 2. Experimental design for transcriptome analysis of cranial neural crest cells and loss-of-function studies. (A) Gastrula stage embryos were 
electroporated with Sox70e2:eGFP construct and incubated until HH1 0. In preparation for FACS, embryos were dissected and enzymatically dissociated. 
GFP+ cells were sorted with 7-AAD exclusion to eliminate damaged or dead cells. RNA from GFP+ and GFP- cells was extracted, reverse transcribed, and 
linearly amplified for high-throughput sequencing. (8) Loss-of-function studies in chicken embryos. Translation-blocking morpholinos (MO) targeted to 
known neural crest regulators were injected into one side of HH6 embryos, with the uninjected side used as an internal control. Electroporated embryos 
were allowed to develop until HH10 and expression of novel neural crest specif ic genes were analyzed by in situ hybridization (ISH). (C) qPCR confirm up- 
regulation of neural crest markers in sorted 5ox70£2:eGFP+ cells in relation to GFP- cells. 



LM04, and nuclear receptor RXRG are expressed only in the neural 
crest at the stages analyzed, in a similar manner to SOX 10 and 
SNAI2. Receptors such as LTK also have sharp and specific ex- 
pression patterns in the neural crest. Transcripts of enzymes like 
BMP2K and matrix molecules such as COL9A3 were expressed 
strongly in the CNC, corroborating our sequencing data. Pre- 
viously, genes up-regulated in the neural crest (Gammill and 
Bronner-Fraser 2002; Adams et al. 2008) were identified by juxta- 
posing neural and non-neural ectoderm to induce neural crest and 
screening macroarrayed embryonic libraries. This process resulted 
in identification of only approximately a dozen transcriptional 
regulators enhanced in the neural crest. In contrast, analysis of 
pure populations of migrating CNC revealed nearly 100 novel 
transcriptional regulators. 

Supplemental Table 4 lists the 50 transcriptional regulators 
with highest enrichment in CNC cells. It includes transcription 
factors, cofactors, modulators, and genes involved in epigenetic 
regulation, 14 of which were previously studied in the context of 
neural crest formation. In total, we estimate that ~ 15%-20% of the 
genes were investigated previously in the context of neural crest or 
craniofacial development or are part of pathways that are known to 
be involved in these processes. Integration of the remainder in the 
CNC GRN will be key for uncovering novel regulatory modules 
that are important for neural crest identity and behavior. 

Effects of perturbation of neural plate border and neural 
crest specifier genes on novel targets 

To identify novel links in the CNC gene regulatory network, we 
assayed the effects of knockdown of known neural plate border 
and neural crest specifier genes on identified targets in the migrating 
CNC population (Fig. 2B). Four genes were chosen as targets for this 



analysis: the nuclear receptor RXRG, the tyrosine receptor LTK, the 
structural protein COL9A3, and the adaptor protein LM04 (Ochoa 
et al. 2012; Ferronha et al. 2013), all of which are robustly expressed 
in the migratory CNC. In addition, we examined effects on tran- 
scription factor EBF1. Seven genes known to be expressed in the 
premigratory neural crest were knocked down using electroporation 
of translation-blocking morpholinos: ETS1, SOX9, PAX7, TFAP2A, 
SOX10, MSX1/2, and MYB. 

Knockdown of neural crest specifiers, ETS1 and 50X9, affected 
the four main genes analyzed (Fig. 6; Supplemental Fig. 1), indicating 
a crucial role for these regulators in CNC development. Indeed, both 
genes have been identified as direct inputs of the SoxlOE2 enhancer 
(Fig. 1; Betancur et al. 2010b). PAX7 is another important gene in the 
neural crest GRN, required for neural crest specification (Basch et al. 
2006) by directly regulating FOXD3 (Simoes-Costa et al. 2012). PAX7 
is maintained in the migratory CNC and coexpressed with SOX10 
(Supplemental Table 4), and its knockdown causes loss of RXRG, LTK, 
and COL9A3 (Fig. 6A), but not of LM04. 

Loss of function of SOX 10 and TFAP2A, key regulators of 
neural crest development (Rada-Iglesias et al. 2012), similarly af- 
fected transcription of LM04 and RXRG but did not alter LTK or 
C0L9A3 expression. Remarkably, loss of the neural plate border 
specifier MSX1/2 resulted in up-regulation of three of the genes 
examined, despite its requirement for neural crest specification 
(Simoes-Costa et al. 2012). The effects of knockdown of neural 
crest specifiers on transcription factor EBF1 (Fig. 2) yielded com- 
pletely different results than the other CNC genes analyzed, sug- 
gesting that its regulation may be indirect (Supplemental Fig. 1). 
MYB knockdown had no significant effect in any of the four 
main CNC transcripts analyzed. Taken together, our results 
point to ETS1 and SOX9 as key regulators of CNC development 
and show differential regulation among the CNC genes surveyed, 
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Figure 3. Cranial neural crest transcriptome data set reveals hundreds of genes up-regulated in cranial neural crest cells. (A) Our transcriptome analysis 
identified 1 301 genes significantly up-regulated in the migratory cranial crest, and 257 genes that were excluded from this population. (6) Genes enriched 
in CNC. The largest fraction are enzymes, followed by transcriptional regulators and transporter molecules. (C) Transcription factors known to be im- 
portant players in the migratory neural crest are among the enriched genes in the CNC transcriptome data set. (D) Most enriched and depleted genes in 
the cranial neural crest. SOX10 and SOX8 are among the genes that are the most enriched in the cranial neural crest, whereas several HOX genes are absent 
from this anterior cell population, as expected. 



highlighting the complex nature of the neural crest gene regulatory 
network. 

Discussion 

Transcriptional regulation of a dynamic cell population like the 
neural crest is an intricate process. The constantly shifting envi- 
ronment to which the migrating CNC cells are exposed, together 
with the onset of their differentiation programs, is likely to con- 
tribute to a high level of regulatory complexity. By greatly 
expanding the candidate genes that can be interrogated for their 
connectivity in the CNC GRN, our data set provides a platform for 



analysis necessary for completion of the migrating CNC GRN. We 
have uncovered approximately 100 novel transcriptional regula- 
tors, many of which undoubtedly play an important role in de- 
termining neural crest identity, potential, and behavior. 

Extensive migratory ability is a defining feature of the neural 
crest. Accordingly, we observe a large number of molecules known 
to play critical roles in cell migration in our data set such as mol- 
ecules involved in changes in cell-cell adhesion and membrane 
receptors. For example, we noted up-regulation of Rho family 
small GTPases (Fig. 4E), important for planar cell polarity and 
mediating cytoskeleton dynamics (Clay and Halloran 2011). In 
fact, nearly half of the up-regulated genes in migrating CNC were 



Genome Research 285 



www.genome.org 



Simoes-Costa et al. 




Fold Change 



Fold Change 



E 



1 Biological Pathway I 


Genes involved H 


Biological Process 1 


Genes involved 


Integrin Signaling 


26 


RNA Processing 


57 


CXCR4 Signaling 


22 


Ubiquitination of Protein 


42 


Protein Llbiquination Pathway 


46 


Migration of Cells 


164 


Signaling by Rho Family GTPases 


33 


Morphology of the Face 


25 


ERK/MAPK Signaling 


34 






Actin Cytoskeleton Signaling 


22 


Disease ■ 


Genes involved 


IGF1 Signaling 


19 


Craniofacial abnormality 


45 


Endothelin 1 Signaling 


27 


Neuroblastoma 


13 


JAK/STAT Pathway 


13 


Osteochondrodysplasia 


10 



Figure 4. Pathways, processes, and diseases associated with the CNC transcriptome data set. Biological pathway analysis shows enrichment of mi- 
grating cranial crest genes involved in integrin (A) and CXCR4 (B) signaling, as expected since both signaling systems are essential for neural crest 
migration. This analysis also reveals potential roles for post-translational regulation in CNC development as a large number of molecules involved in RNA 
splicing and processing are up-regulated (C). (D) Genes involved in diseases such as osteochondrodysplasia are present in the data set, which is consistent 
with the contribution of the neural crest to the craniofacial skeleton. (£) Enrichment of CNC genes in other pathways, processes, and diseases. 



enzymes. This is likely due to the constant structural rearrange- 
ments, extensive proliferation, and active signal transduction, 
such as the MAPK/ERK pathway, ongoing in migratory CNC cells 
(Fig. 4; Newbern et al. 2008; Stuhlmiller and Garcia-Castro 2012). 
Furthermore, numerous enzymatic epigenetic modifiers (Hu et al. 
2012; Strobl-Mazzulla and Bronner 2012) including helicases and 
topoisomerases are active in the CNC. 

Consistent with classical experiments revealing differences in 
developmental potential between cranial and trunk neural crest, 
our data set shows enrichment of major components of the 
chondroblast differentiation program in migrating CNCs. These 
include SOXD family members, MIA/CD-RAP, TFAP2A transcrip- 
tion factors, and cartilage-specific matrix molecules COL2A1, 
COL9A3, COL11A1 and ACAN (Supplemental Table 1). The pres- 
ence of these transcripts may reflect early specification of a subset 



of CNC toward an ectomesenchymal fate, perhaps explaining why 
some cranial crest cells differentiate into ectopic cartilage after 
grafting to the trunk (Nakamura and Ayer-le Lievre 1982; Le 
Douarin et al. 2004). Comparative genomic analysis between dis- 
tinct neural crest subpopulations holds the promise of revealing 
key circuits that endow the CNC with its unique features. 

The complexity of the CNC regulatory program is highlighted 
in our functional analysis, which reveals many previously un- 
known connections between upstream regulators and the identi- 
fied effector genes (Fig. 7). Results reveal important trends in reg- 
ulatory control of CNC-specific transcripts. For instance, our data 
show that ETS1 and SOX9 are required for the expression of four of 
five CNC genes analyzed. This suggests a model in which ETS1 and 
SOX9 may act as pan-neural crest regulators during CNC migration 
(Fig. 7), hinting at a common subcircuit controlling expression of 
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Figure S. In situ hybridization screen and qPCR tests validate CNC transcriptome data set. In situ hybridization confirms enrichment of migrating CNC 
genes. These include well-known neural crest markers such as 5OX70, ETS1, and SNAI2 (A-Q as well as novel genes previously not known to be expressed 
in the migrating CNC (D-T). (U) qPCR tests corroborate enrichment found in the RNA-seq data. Vertical error bars on data points represent the standard 
error of mean obtained from replicates. 



CNC-specific transcripts. This view is supported by ris-regulatory 
analysis of bona-fide neural crest markers (Betancur et al. 2010b; 
Simoes-Costa et al. 2012; Barembaum and Bronner 2013), identi- 
fying SOX9 and ETS1 as direct inputs in neural crest enhancers. 

Apart from common inputs from SOX9 and ETS1, each of the 
target genes analyzed requires a different combination of upstream 
regulators for transcriptional activation. This is not unexpected 
since the CNC is complex and its regulation involves the interplay 
of a large number of transcription factors (Betancur et al. 2010a). 
Since the four genes analyzed have very similar expression pat- 



terns, our results suggest the existence of concurrent genetic pro- 
grams within the CNC that result in similar transcriptional out- 
comes. According to this model (Fig. 7), different combinations of 
transcription factors would result in transcriptional activation of 
distinct genes in CNCs. SOX9 and ETS1 would thus be a part of a 
common essential genetic circuit working in conjunction with dif- 
ferent regulators (PAX7, SOX10, TFAP2A) to activate transcription. 

Interestingly, the same regulator can have distinct effects at 
different developmental times and in different modules of the 
CNC GRN. For instance, MSX1 plays a critical early role as an ac- 
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Figure 6. Regulatory interactions between classic and novel regulators of neural crest development. Functional studies reveal epistatic interactions 
between neural plate border and neural crest specif iers with CNC genes. (A) Summary of the results obtained from the loss-of-function experiments. (B-F) 
Representative embryos for each of the phenotype classes: (6) no phenotype; (C) mild knockdown; (D) strong knockdown; (£) mild up-regulation; and (F) 
strong up-regulation. (C) Quantification of knockdown experiments. 



tivator of neural crest specification (Monsoro-Burq et al. 2005; 
Simoes-Costa et al. 2012), but our data show that in the migrating 
CNC, it modulates or represses gene expression. Consistent with 
repressor function, MSX1 has been shown to recruit Polycomb 
complexes during embryonic development to shut down target 
genes during myoblast differentiation (Wang et al. 2011). A likely 
explanation is that the same regulator often can be used re- 
iteratively but to different effect depending upon available partners. 
These results highlight the complexity of regulatory interactions in 
the neural crest GRN and demonstrate that the consecutive pro- 
grams that drive neural crest gene expression can be quite distinct, 
even if they utilize some of the same suite of regulators. 

The migrating CNC transcriptome data set presented here pro- 
vides a platform for a thorough exploration of the transcriptional state 



of CNC cells during their active phase of migration. The data set will 
facilitate expansion of the terminal modules of the neural crest GRN. 
This will lead to a better understanding of how the genetic program 
controlling neural crest identity unfolds and is converted into discrete 
regulatory modules that drive differentiation into a variety cell 
types. Our results point to a variety of molecules and pathways that 
have yet to be explored in the context of neural crest development 
and offer a panoramic view of the CNC transcriptional landscape. 

Methods 

Electroporation of enhancers and time-lapse imaging 

Chicken embryos were incubated to Hamburger and Hamilton 
(HH)4 (Sauka-Spengler and Barembaum 2008). Embryos were dis- 
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Figure 7. Migrating cranial crest data set allows for expansion of neural crest gene regulatory network and identification of new regulatory links. 
Transcriptome analysis of the migrating cranial neural crest expands the number of regulators in the migratory module of the CNC gene regulatory 
network. Loss-of-function experiments identify new links between well-characterized regulators and novel neural crest effectors. 



sected (Chapman et al. 2001), placed in Ringer's solution, injected 
with 2 ug/uL of the pTK-Soxl0e2:eGFP construct (Betancur et al. 
2010b), electroporated, and placed in individual culture dishes 
containing albumen supplemented with Penicillin-Streptomycin 
(Sigma). For dynamic analysis of the SoxlOE2 enhancer, HH4 chick 
embryos were electroporated as described above. After ~8 h of in- 
cubation, embryos were imaged as described (Ezin and Fraser 2008). 

Tissue dissociation and FACS 

After electroporation, embryos were incubated at 37°C until 
HH10 (—14-16 h of incubation) and dissected above the level of 
rhombomere 3 to avoid contamination of otic GFP+ cells. Dis- 
sected heads were dissociated in dispase/trypsin. Reaction was 
stopped by addition of 10 times the volume of Hank's solution. 
Cells were washed, passed through a 40u.m strainer (BD) and 
resuspended in Hank's solution plus 0.5% BSA, 5 units/mL DNase I 
(Roche). Cells were sorted using BD FACSAria Cell Sorter (BD Bio- 
sciences) with 7-AAD exclusion to eliminate dead/damaged cells. 
GFP+ cells from dissected heads and GFP- cells from whole em- 
bryos (WE) were used for library construction and sequencing. 

RNA Extraction and Amplification 

Sorted cells were pelleted, washed with PBS, and immediately lysed 
with Lysis Buffer. RNA was extracted with RNAqueous-Micro Kit 
(Ambion AM 1931); genomic DNA was removed by 30 min of 
TURBO DNase treatment. Quality of the RNA was assayed in an 
Agilent 2100 Bioanalyzer. All samples had RIN>8. RNA was sub- 
sequently amplified with the Ovation RNA-Seq System V2, yield- 
ing double-stranded cDNA that was used for library construction. 

Library building and sequencing 

Sequencing libraries were built according to Illumina Standard 
Protocols. SR50 sequencing was performed in a HiSeq Illumina 
machine. Sequence reads were aligned to the Gallusgallus genome 



with Bowtie (version 0.12.7) and TopHat (version 1.3.1.0SX_ 
x86_64). The following reads and percentage of reads maps were 
obtained for each of four libraries: CNC1 (21.5 million reads/87% 
mapped), CNC2 (26.2 million reads/50% mapped), WEI (21.7 
million reads/81% mapped) and WE2 (30.7 million reads/64% 
mapped). Cufflinks and Cuffdiff were used to calculate gene expres- 
sion levels and identify differentially expressed transcripts. Two CNC 
and two whole embryo samples were independently collected, se- 
quenced and mapped, and subsequently analyzed in Cufflinks 
and Cuffdiff as replicates. The Pearson correlation coefficients 
obtained between replicates were 0.925 for the CNC samples and 
0.957 for the WE samples. A gene was considered to be expressed by 
the CNC transcriptomic data set if (1) the RPKM values (Reads Per 
Kilobase of exon per Million fragments mapped) for that gene were 
above 5 RPKM in the population with higher expression levels; 
(2) difference in expression levels was significant according to sta- 
tistical analysis; and (3) the gene had at least twofold up-regulation/ 
depletion in the CNC. The expression levels on the CNC data set 
(Supplemental Table 1) contain values for samples that were 
processed as replicates in TopHat. Biological pathway analysis was 
performed with the aid of Ingenuity iReports, after conversion of the 
chicken Ensembl gene models to mouse orthologous counterparts. 

cDNA synthesis and qPCR 

RNA from GFP+ and GFP- sorted cells was isolated using the 
RNAqueous-Micro Kit (Ambion AM1931), according to the manu- 
facturer's protocol. cDNA was synthesized using Superscript III First 
Strand Synthesis kit (Invitrogen) with random primers for RT prim- 
ing. qPCR was performed using SYBR Green (Bio-Rad) according to 
the manufacturer's instructions. For reference genes, we used 
qPCR primers designed for the HPRT1 and TFRC genes. 

In situ hybridization and antibody labeling 

Whole mount in situ hybridization was performed as described 
(Acloque et al. 2008). Embryos electroporated with the SoxlOE2: 
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eGFP construct were labeled using HNK-1 and SOX10 antibodies 
(1:100), and goat anti-mouse IgM Alexa 350 or anti-mouse IgG 
Alexa 350, respectively (1:1000; Molecular Probes). Whole mount 
immunostaining was performed as described (Ezin et al. 2009). For 
in situ hybridization, fragments of candidate genes ranging from 
600-800 were PCR cloned using the TOPO TA Cloning Kit from 
Invitrogen (K4610-20). Clones were sequenced and used as tem- 
plates for DIG-labeled RNA probes. Expression of approximately 80 
candidate genes in the CNC was verified in chicken embryos 
(Supplemental Table 1). 

Loss-of-function studies 

For morpholino knockdown, unilateral electroporation was per- 
formed in HH6 chicken embryos as described (Barembaum et al. 
2000), using the following morpholinos: Msxl, Pax7, and Etsl 
(Simoes-Costa et al. 2012); c-Myb, Sox9, and SoxlO (Betancur et al. 
2010b); andTafp2a (Barembaum and Bronner 2013). Embryos were 
collected and fixed in 4% paraformaldehyde for 1 h and screened 
under fluorescence for efficient morpholino incorporation. 

Data access 

All sequencing data have been deposited in the NCBI BioProject 
database, under BioProject PRJNA226749 (http://www.ncbi.nlm. 
nih . go v/biopro j ect) . 
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